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ABSTRACT 



Context. The elfect of a frequency dependent shift of the VLBI core position (known as the "core shift") was predicted more than 
three decades ago and has since been observed in a few sources, but often within a narrow frequency range. This effect has important 
astrophysical and astrometric applications. 

Aims. To achieve a broader understanding of the core shift effect and the physics behind it, we conducted a dedicated survey with 
NRAO's Very Long Baseline Array (VLBA). 

Methods. We used the VLBA to image 20 pre-selected sources simultaneously at nine frequencies in the 1.4—15.4 GHz range. The 
core position at each frequency was measured by referencing it to a bright, optically thin feature in the jet. 

Results. A significant core shift has been successfully measured in each of the twenty sources observed. The median value of the 
core shift is found to be 1.21 mas if measured between 1.4 and 15.4 GHz, and 0.24 mas between 5.0 and 15.4 GHz. The core position, 
r^, as a function of frequency, y, is found to be consistent with an oc y"' law. This behavior is predicted by the Blandford & Konigl 
model of a purely synchrotron self-absorbed conical jet in equipartition. No systematic deviation from unity of the power law index 
in the r^iv) relation has been convincingly detected. 

Conclusions. We conclude that neither free-free absorption nor gradients in pressure and/or density in the jet itself and in the 
ambient medium surrounding the jet play a significant role in the sources observed within the 1.4-15.4 GHz frequency range. These 
results support the interpretation of the parsec-scale core as a continuous Blandford- Konigl type jet with smooth gradients of physical 
properties along it. 

Key words. Galaxies: active - galaxies: jets - quasars: general - radio continuum: general 



1. Introduction 

In VLBI images of relativistic AGN jets, the "core" is the 
most compact feature near the apparent base of the jet. It is 
identified with the surface at which the optical depth t,, ^ 1 
(photosphere) in a contin uous jet flow wi th gradually changing 
physical properties (iBlandford & Konigll 1 1979b . The position 
of this surface is a function of the observation frequency, v, 
and the exact form of this function depends on the absorption 
mechanism. The above interpretation of the core is supported by 
observations of the "core shift" effect serendipitously discovered 
bv lMarca ide & Shapir d(ll984l) and l ater in ve stigated by , amon g 



others, Zensus Cohen. & UnwinI (19951). LobanoyI ( 



iKovalev et alJ (l2008h . and lO' Sullivan & Gabuzdal (120091) . 



1998), 



Additional support comes from observations of the related 
increase in the apparent core size toward lower frequencies 
(lUnwin et al.lll994lYang et~ani2008h . 

An alternative interpretation of the VLBI core as a stand- 
ing shock in the jet is discussed by Marscher (2006, 2008). This 
interpretation may be relevant for high-frequency observations 
of some sources. The standing shock model does not predict a 
shift in the apparent core position with frequency unless one re- 
sorts to an assumption that there is a series of multiple standing 
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shocks, and different shocks are observed as the core at different 
frequencies. In some objects the core could be a region where 
the jet bends such that it is more aligned with the line of sight 
leading to an increased Doppler beaming factor relative to jet 
regions farther upstream (lMarscherll2006l) . It has even been dis- 
cussed sometimes, that the radio core may not be part of the jet 
at all. 

The most likely absorption mech anism acting in the com- 
pact jet is synchrotron self-absoiption (lKoniglll981l) . If the jet is 
freely expanding and there is an equipartition between the parti- 
cle and magnetic field energy densities, the position (v) of the 
Ty = 1 surface will be rdv) oc y^'. If the above assumptions do 
not hold (e.g., synchrotron self-absorption is not the dominating 
absorption mechanism), we can expect a different dependency 
of rdv). 

Investigation of the core shift effect is important for gaining 
deeper understanding of the structure and physical conditions in 
ultracompact AGN jets. It may also provide information about 
the pressure and density gradients in ambient medium surround- 
ing VLBI-scale jets - the broad-line region ( BLR) and the in- 
ner part of the naiTow-line region (NLR) (e.g., Eobanov 2007!). 
The impact of the effect on interpretations of current and fu- 
ture radio Very Long Baseline Interfero metry (VLBI) astro met- 
ric measureme nts has been discussed bv lRioia et al. I (l2005h and 
iPorcasI (l2009h . The core shift is expected to introduce a sys- 
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Table 1. The large core shift source sample observed with the VLBA 



lAU name Alias R.A. (J2000) Dec. (J2000) z Optical class VLBA epoch 



0148+274 




01:51:27.146174 


+27:44:41.79363 


1.26 


QSO 


2007- 


-03-01 


0342+147 




03:45:06.416545 


+ 14:53:49.55818 


1.556 


QSO 


2007- 


■06-01 


0425+048 


OF 42 


04:27:47.570531 


+04:57:08.32555 


0.517" 


AGN 


2007- 


■04-30 


0507+179 




05:10:02.369133 


+ 18:00:41.58160 


0.416 


AGN 


2007- 


■05-03 


0610+260 


3C 154 


06:13:50.139161 


+26:04:36.71971 


0.580 


QSO 


2007- 


■03-01 


0839+187 




08:42:05.094175 


+ 18:35:40.99050 


1.272 


QSO 


2007- 


■06-01 


0952+179 




09:54:56.823616 


+ 17:43:31.22204 


1.478 


QSO 


2007- 


■04-30 


1004+141 




10:07:41.498089 


+ 13:56:29.60070 


2.707 


QSO 


2007- 


■05-03 


1011+250 




10:13:53.428771 


+24:49:16.44062 


1.636 


QSO 


2007- 


■03-01 


1049+215 




10:51:48.789077 


+21:19:52.31374 


1.300 


QSO 


2007- 


■06-01 


1219+285 


W Comae 


12:21:31.690524 


+28:13:58.50011 


0.161'' 


BLLac 


2007- 


■04-30 


1406-076 




14:08:56.481198 


-07:52:26.66661 


1.493 


QSO 


2007- 


■05-03 


1458+718 


3C 309.1 


14:59:07.583927 


+71:40:19.86646 


0.904 


QSO 


2007- 


■03-01 


1642+690 




16:42:07.848505 


+68:56:39.75639 


0.751 


QSO 


2007- 


■04-30 


1655+077 




16:58:09.011464 


+07:41:27.54034 


0.621 


QSO 


2007- 


■06-01 


1803+784 




18:00:45.683905 


+78:28:04.01839 


0.680 


QSO 


2007- 


■05-03 


1830+285 




18:32:50.185622 


+28:33:35.95514 


0.594 


QSO 


2007- 


■03-01 


1845+797 


3C 390.3 


18:42:08.989895 


+79:46:17.12825 


0.056 


AGN 


2007- 


■06-01 


2201+315 




22:03:14.975788 


+31:45:38.26990 


0.298 


QSO 


2007- 


■04-30 


2320+506 




23:22:25.982173 


+50:57:51.96364 


1.279 


QSO 


2007- 


■05-03 



Column designation: Col. 1 - lAU source name (B1950), Col. 2 - commonly use d source name, Cols. 3 and 4 - VLBI position, for de- 
tails see http://astrogeo.org/vlbi/solutions/rfc_281Qc/ and lBeaslev et al.l ll2b02'). Fomalont et al. (2003), Petrov et al. (2005, 200(1), 
iKovalev et ali l2007i) . iPetrov et aU ^2008, 2009), Cols. 5 and 6 - redshift and optical class from iVeron-Cettv & Veronl(i201Q) . Col. 7 ~ epoch of 
multifrequency VLBA observations reported in this paper 
" Spectroscopic redshift obtained by Afanasiev et al .1(120031) . 
Photometric redshift, see finke et al.. (.2008.) . 



tematic offsets of approximately 0.1 mas between the optical 
(6000 A) and radi o (8 GH z) positions of reference extragalac- 
tic sources ( Koval ev et al ]!2008), which needs to be taken into 
account for accurate radio-optical reference frame al ignment in 
the er a of modern spac e-based astrometric m issions (iLindegrenl 
l2007h such as GAIA (iPerrvman etani200lh and for accurate 
spacecraft navigation with VLBI (e.g.. [ Hildebrand et alJll994 
[Sekido et al. 2004, Pogrebe nkoet al.l |2009|). 

In this paper, we present observational results from a dedi- 
cated nine-frequency (1.4-15.4 GHz) VLBA survey of the core 
shift effect in 20 compact extragalactic radio sources. The se- 
lection of targets, the VLBA observation setup and our data re- 
duction strategy are described in Section |2] Section |3] presents 
the core shift measurement technique. In Section |4] we discuss 
the results and their implication for the dominating absorption 
mechanism in the core region. Section |5] presents a brief sum- 
mary of this work. Throughout this paper, we adopt a A-CDM 
cosmology, with the following values for the co smological pa- 
rame ters: //() = 0.71, flm = 0.27, and = 0.73 (iKomatsu et alj 
120091) . 

2. VLBA observations 

2.1. Sample selection 

A blind survey for sources with measurable core shifts would 
be an extraordinarily time consuming project because it requires 
performing simultaneous multifrequency, high-resolution, high 
dynamic range VLBI observations. Therefore, we turned to an 
existing large-scale geodeti c VLBI datab a se to fi nd promising 
candidates for the survey. iKovalev et al.l (l2008h have imaged 
and analyzed 277 sources from the Research and Development 
VLBA program (RDV, see Fey & Chariot 1997 and Petrov et aL 
12009 ) observations made in 2002-2003. RDV experiments in- 
volve simultaneous dual frequency 2.3/8.6 GHz {S /X band) ob- 



servations with the VLBA and up to ten other globally dis- 
tributed antennas. The data processing technique and imaging 
results are described by Pushkarev & Kovalev (2011). The core 
shift was measured by Kov alev et al.i (2008) by referencing the 
core position to optically thin jet features, the positions of which 
are not expected to change with frequency. 

During the preparation of the IKovalev et al.l ('2008') work, 1 8 
sources with large core shifts between 2.3 and 8.6 GHz were se- 
lected for detailed multifrequency study with the VLBA, which 
we report here. These sources were selected in such a way that 
they not only show large core shifts in this frequency range, 
but are also particularly suitable for measuring it by referencing 
the core position to bright jet features - the technique used in 
this paper. This procedure is needed because the absolute posi- 
tion information is lost during the phase self-calibration process 
necessary for VLBI imaging; therefore, it is not known a priori 
how images at different frequencies should be aligned. Two more 
sources, 2201+315 and 3C 309.1 (1458+718), were ad ded to the 
list in order to continue our previous core shift studies (iLobanovl 
ll998llRos & Lobanovl200"2LlHoman & Kovalevl2003l) . 

2.2. VLBA observation setup 

Twenty sources selected as promising candidates for a detailed 
core shift investigation were observed with the VLBA (iNapieil 
11994 ) during four 24 h sessions in March - June 2007. The list 
of the observed sources is presented in Table [1] Each source 
was observed simultaneously using L, S , C, X, and K,, band re- 
ceivers (according to the Institute of Electrical and Electronics 
Engineers or IEEE nomenclature, see Table |2) in the 1.4- 
15.4 GHz range. In each band four 8 MHz-wide frequency chan- 
nels (IFs, Table ^ were recorded in both right and left circular 
polarizations with 2-bit sampling and a total aggregate bit rate of 
256 Mbits/sec. The correlation of the data was performed at the 
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Table 2. IF central frequencies 



IEEE Band IF Frequency 
(MHz) 



L 


1 


1404.49 


L 


2 


1412.49 


L 


3 


1658.49 


L 


4 


1666.49 


S 


1 


2275.49 


S 


2 


2283.49 


s 


3 


2383.49 


s 


4 


2391.49 


c 


1 


4604.49 


c 


2 


4612.49 


c 


3 


4999.49 


c 


4 


5007.49 


X 


1 


8104.49 


X 


2 


8112.49 


X 


3 


8425.49 


X 


4 


8433.49 


K„ 


1 


15353.49 


K„ 


2 


15361.49 


K„ 


3 


15369.49 


K„ 


4 


15377.49 



Column designation: Col. 1 - Radio band name according to the IEEE 
radar band nomenclature, Col. 2 - number of the frequency channel 
(IF), Col. 3 - central frequency of the frequency channel (IF). 



VLBA Array Operation Center in Socorro, NM, USA, with an 
averaging time of 2 seconds. Data in L, 5, C, and X bands were 
divided into two sub-bands (two IFs in each sub-band) centered 
at 1.4, 1.7, 2.3, 2.4, 4.6, 5.0, 8.1, and 8.4 GHz, respectively. The 
K,i band was not divided into sub-bands in order to achieve sim- 
ilar image sensitivity to other frequencies in a comparable inte- 
gration time. All four IFs in the K,, band were stacked together 
around 15.4 GHz. The central frequencies of the sub-bands were 
chosen in such a way that at least one sub-band in each band is 
centered at a frequency for which antenna gain curve measure- 
ments are available. A special procedure was applied to ensure 
accurate amplitude calibration of all sub-bands (Section l23T l. 

2.3. VLBA data calibration and imaging 

The initial calibration was conducted in AIPS dGreise ril l2003h 
following the standard VLBA calibration procedure involving 
(i)a priori amplitude calibration with measured antenna gain 
curves and system temperatures, ( // jphase calibration using the 
phase-cal signal injected during observations, and (iii)fnnge fit- 
ting. The fringe fitting was performed by the task FRING. A 
separate solution for the group delay and phase rate was made 
for each frequency channel (IF). As the final step in calibration, 
bandpass corrections were applied utilizing the task BPASS. 

The sources were imaged independently in each band, us- 
ing the CLEAN al gorithm ( Hogbom 1974) implemented in the 
Difmap software ( Shepherdll997 ). Global amplitude corrections 
for left and right circular polarization for each IF at each an- 
tenna were determined by comparing the total intensity CLEAN 
model obtained with the initially calibrated data (again using 
Difmap). The amplitude corrections obtained were averaged 
over all the sources observed in the experiment with the excep- 



Table 3. Amplitude corrections for the BK134 VLBA experi- 
ment 



Ant. 


Band 


IF 


Epoch 


Polarization 


Corr 


BR 


L 


1 


All 


LCP 


1.17 


BR 


L 


2 


All 


LCP 


1.13 


LA 


L 


1 


All 


LCP 


0.89 


LA 


L 


2 


All 


LCP 


0.90 


HN 


L 


1-4 


2007-06-01 


LCP RCP 


0.75 


OV 


L 


1 


All 


RCP 


1.17 


OV 


L 


2 


All 


RCP 


1.15 


BR 


S 


3 


All 


LCP 


1.13 


LA 


S 


1 


All 


LCP 


1.12 


LA 


S 


2 


All 


LCP 


1.09 


LA 


S 


3 


All 


LCP 


1.19 


LA 


S 


4 


All 


LCP 


1.23 


HN 


S 


3 


All 


LCP 


0.90 


HN 


S 


4 


All 


LCP 


0.89 


KP 


S 


1 


All 


LCP 


1.09 


KP 


S 


2 


All 


LCP 


1.13 


NL 


s 


3 


All 


LCP 


1.09 


MK 


S 


3 


All 


LCP 


1.13 


MK 


S 


4 


All 


LCP 


1.14 


LA 


S 


1 


All 


RCP 


0.85 


LA 


S 


2 


All 


RCP 


0.85 


LA 


S 


3 


All 


RCP 


0.76 


LA 


S 


4 


All 


RCP 


0.79 


OV 


S 


3 


All 


RCP 


0.86 


OV 


s 


4 


All 


RCP 


0.89 


KP 


c 


1 


All 


LCP 


1.10 


KP 


c 


2 


All 


LCP 


1.11 


BR 


c 


1 


All 


LCP 


0.91 


KP 


c 


1 


All 


RCP 


1.11 


KP 


c 


2 


All 


RCP 


1.09 


MK 


c 


1 


All 


RCP 


1.24 


MK 


c 


2 


All 


RCP 


1.17 


OV 


X 


1 


All 


LCP 


1.21 


OV 


X 


2 


All 


LCP 


1.20 


KP 


X 


1 


All 


RCP 


0.91 


KP 


X 


2 


All 


RCP 


0.93 


BR 


K„ 


1 


All 


LCP 


0.94 


BR 


Ku 


4 


All 


LCP 


0.95 


SC 


K„ 


1 


All 


LCP 


1.05 


KP 


Ku 


2 


All 


RCP 


1.42 



Column designation: Col. 1 - Antenna name, Col. 2 - radio band 
name, Col. 3 - number of the frequency channel (IF), Col. 4 - epoch to 
which the amplitude correction was applied, Col. 5 - polarization (right 
or left circular). Col. 6 - amplitude correction coefficient. 



tion of 0610-H 260and 1830+ 285. These two sources have steep 
radio spectra (iKovalev et al.lil999, .2002) and were not used for 
the amplitude corrections computation in order to avoid intro- 
ducing systematic errors into the corrections. The significant am- 
plitude correction factors that were found are listed in Table [3] 
These amplitude corrections were introduced into the dataset us- 
ing the AIPS task CLCOR. After that L,S,C, and X band data 
were split into two sub-bands as described above, and final to- 
tal intensity imaging for each sub-band was conducted indepen- 
dently using Difinap. Naturally weighted CLEAN images of the 
observed sources between 1.4 and 15.4 GHz are presented in 
Fig. [U available in the electronic version of the paper. 
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The typical VLBA amplitude calibration accuracy in the 
1-10 GHz frequency range is estimated to be ~ 5 % accord- 
ing to the VLBA Observational Status Summar^Q. Kovalev et al. 
(I2005h and lPetrov et alJ (l2008h have confirmed this estimate to 
be correct for S, X, and K,, bands by comparing VLBA in- 
tegrated flux density with a single-dish flux density measured 
from 1 to 22 GHz quasi-simultaneously at RATAN-600 (see 
description of the program in Kovalev et al. 1999, 2002) and 
the data from the University of Michigan Radio Astronomy 
Observatory (UMRAO ) moni toring program (lAUer et al.lll985[ 
lAller. Allen & HughesI 1 19921 |2003). This agrees wifli flie 
VLBA Ku band a mplitude calibration accuracy reported by 
iHoman et al.l (12002 ). We checked and confirmed the above es- 
timates for this particular dataset by using our RATAN-600 
data at 1-22 GHz and VLA measurements at 5 and 8 GHz 
from the VLA/VLBA Polarization Calibration PrograrcQ- This 
check was particularly efficient with our slowly varying calibra- 
tor 09234-392 observed in each of the four VLBA sessions. It has 
negligible kiloparsec-scale emission, which makes this compar- 
ison possible. 



3. Model fitting and core shift measurements 

The structure of each source at each frequency was modeled 
in the visibility (mv) plane with a number of circular Gaussian 
components using the Difinap software (iShepherdI 1 997h . A suit- 
able number of components (typically ~ 6) was used to account 
for all significantly detected features of the source structure at a 
given frequency. The components were visually cross-identified 
between frequencies on the basis of their position with respect 
to the core. The VLBI core is identified at each frequency as 
the bright component at the apparent jet base. As pointed out in 
Section ITTl the target sources were preselected in a way that 
such an identification would be possible in a wide frequency 
range for the brightest jet components. 

One bright, optically-thin jet component, which could be 
identified across all (or most of) the observing frequency range, 
was chosen, and its distance from the core was measured at each 
frequency. Model component positions and sizes for the core, 
as well as the jet feature used for this analysis, are indicated in 
Fig. [T] In some rare cases we were unable to construct a con- 
sistent model to allow for accurate and robust component cross- 
identification with other frequencies. In these cases no model 
components are shown in Fig. [T] and no data are used in the 
analysis. Figure |2] shows the distribution of spectral indexes of 
the reference components, indicating that they are optically thin 
in the considered frequency range. 

The technique of using an optically thin jet component as 
a reference point for the core position measurement has pre- 
viously been s uccess fully tested on S /X band RDV data by 
iKovalev et alj j2008l) in several ways: (i) similar core shift 
values were obtained referencing to different jet features of 
the same source; (ii) the core shift for the quasar 1655-1-077 
agreed wit hin the errors with previo us phase-referenced mea- 
surements (Homan & Kovalevl l2003l) : (Hi) two RDV epochs of 
the q uasar 164 2-f690 resulted in values being in good agree- 
ment. iKovalev et all (f2008l) have also investigated the effect of 
the different wv-coverage resulting in a different resolution (and 
amount of blending between the core and nearby jet regions) at S 
and X bands on core shift measurements. The effect was found 



' http : //www . vlba . nrao . edu/astro/obstatus/current/] 
^ jhttp : //www ■ aoc . nrao . edu/~ smyers/calibration/| 
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Fig. 2. Distribution of the spectral index (a, flux density 5,, oc v") 
of reference jet components. The indexes were computed by fit- 
ting a power law to the 1.4-15.4 GHz spectrum. The compo- 
nent flux at each frequency was obtained from the wv-model. 
The source with the steepest reference component spectrum 
(0- = -1.29) is 1011-H250. 

to be small compa red to the core posit ion shift between these 
bands (see Fig. 5 in lKovalev et alJl2008l) . 

In addition to that, we performed the following test to look 
for indications of blending in our 1.4-15.4 GHz results. If sig- 
nificant blending is present, we should measure systematically 
larger core shifts for sources with jet directions aligned with the 
position angle (PA.) of the elongated VLBA beam. However, no 
such dependence is found in the data. 

To assess the quality of image alignment using a reference 
jet component, we inspected spectral index maps of the sources. 
The spectral index maps were computed by constructing a data 
cube from the images at different frequencies convolved with the 
same beam (naturally weighted beam corresponding to the low- 
est frequency) and fitting a power law to the spectrum in each 
individual pixel of the data cube. The power law index (corre- 
sponding to the slope of a straight line if plotted in logarithmic 
scale) is the spectral index. A map with a smooth spectral index 
gradient along the jet and an absent/weak gradient across the 
jet would indicate, in general, the correct alignment of images 
obtained at different frequ encies (see, e.g.. lMarcaide & Shapirol 
1 1984 IKovalev et alJl2008h . In all cases the spectral index maps 
revealed no major problems with image alignment, and the core 
was found to be the only optically thick feature in the jet. The 
detailed discussion of the spectral imaging results will be pre- 
sented elsewhere. 

An alternative approach to multifrequency VLBI image 
alignment ba sed on two-di r nensi onal cross-correlation has been 
proposed by IWalker et~an (|2000). It ha s been successfuUy ap- 
plied for core shift measurements by lO' Sullivan & Gabuzdal 
d2009h in sources with smooth jets where referencing to a dis- 
tinct bright jet component was not possible. A potential disad- 
vantage of this method is that, in principle, it may introduce ar- 
tificial offsets between frequencies in the presence of spectral 
gradients along the jet. A similar problem may affect the model 
component-based method used here, if there is a large spec- 
tral gradient across a resolved reference component. However, 
the spectral index gradient is smaller for an individual compo- 
nent then for a large section of parsec-scale jet used in cross- 
correlation analysis. Early observations of the frequency depen- 
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Fig. 3. Separation of the core from a reference optically-thin jet component as a function of frequency. The curve represents the 
best-fit function rdv) - a + bjv. The best-fit parameters are presented in Table |5] 
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Fig. 3. continued... 

Table 4. Core shift in milliarcseconds with respect to its position measured at 15.37 GHz 











Frequency (GHz) 








1.41 


1.66 


2.28 


2.39 


4.61 


5.00 


8.11 


8.43 


Name 


















0148+274 


1.41 


1.30 


1.10 


1.16 


0.34 


0.37 


0.14 


0.10 


0342+147 


0.71 


0.64 


0.30 


0.28 


0.07 


0.07 


0.03 


0.06 


0425+048 


1.36 


1.11 


0.98 


0.82 


0.20 


0.12 


0.17 


0.10 


0507+179 


0.90 


0.98 


0.60 


0.65 


0.03 


0.08 


0.00 


0.00 


0610+260 


2.20 


1.99 


1.53 


1.58 


0.45 


0.50 


0.21 


0.22 


0839+187 


2.46 


2.15 


1.03 


0.98 


0.32 


0.29 


0.20 


0.22 


0952+179 


1.07 


1.00 


0.59 


0.60 


0.47 


0.44 


0.16 


0.17 


1004+141 


1.61 


1.21 




1.03 


0.31 


0.31 


0.08 


0.05 


1011+250 


1.27 


1.22 








0.20 


0.12 




1049+215 


1.21 


0.91 




0.70 


0.35 


0.28 


0.08 


0.14 


1219+285 




1.18 


0.88 


0.81 


0.21 


0.21 


0.04 


0.00 


1406-076 


0.78 


0.70 


0.48 


0.49 


0.23 


0.22 


0.07 


0.09 


1458+718 


1.39 


1.25 


0.98 


0.94 


0.18 




0.08 


0.08 


1642+690 


1.01 


0.95 


0.37 


0.32 


0.06 


0.04 


-0.06 


-0.06 


1655+077 


0.88 


0.66 






0.14 


0.12 


-0.01 


0.02 


1803+784 


0.60 


0.49 


0.48 


0.46 


0.04 


0.04 


0.09 


0.01 


1830+285 


2.11 




0.84 


0.93 


0.35 


0.37 


0.01 


0.04 


1845+797 




1.49 


0.48 


0.58 


0.22 


0.29 


0.08 


-0.00 


2201+315 






1.17 


1.16 


0.39 


0.42 


0.14 


0.11 


2320+506 


0.62 


0.60 


0.40 


0.45 


0.09 


0.08 


-0.02 


-0.03 


Adopted 


















uncertainty 


0.07 


0.07 


0.07 


0.07 


0.04 


0.04 


0.04 


0.04 



Designation: Col. 1 - lAU source name (B1950), Cols. 2 to 9- core shift measured at the given frequency (in milliarcseconds). The last row 
of the table presents the adopted measurement uncertainty at the given frequency (in milliarcseconds). The adopted measurement uncertainty at 
15.37 GHz is 0.04 mas. 



dent difference in distance between the core and jet components 
were actually interpreted by Biretta, Moore, & Cohen (1986) as 
a sign of spectral gradients across the components. This inter- 
pretation is disfavored by the detection of cor e position shift 
with frequency in phase - referenced experiments ([Marcaide et al] 
I1985L iRioia et all [19971 |Homan & Kovalev' 2003) and the fact 
that strong spectral gradients across individual jet components 
are not routinely observed if images are aligned correctly. 

A significant core shift has been detected in each of the 
twenty sources observed within our VLBA program (Table 0. 
The effect is especially pronounced at frequencies below 5 GHz. 
The results are discussed in detail in Section |4] and presented 
visually in Fig. [3] 

To statistically estimate typical uncertainty cr^ of the mea- 
sured distance r between the core and reference jet component, 
we used pairs of sub-bands of X, C, S , and L bands, assuming 
that the source structure at these closely separated frequencies 
is essentially the same and the non-zero difference 6r between 
sub-bans in the model fit parameters is determined purely by er- 
rors (Fig. |4]i. This assumption turns out to be correct for X, C, 



and S bands, while at L band a significant core shift between 
the sub-bands becomes evident. We underline that this approach 
is possible since data processing of these sub-bands was done 
completely independently. It should be also noted that the core- 
shift effect is found to be not significant between the two sub- 
bands at S -band most probably due to the relatively narrow fre- 
quency range (Table|2]l achievable at this particular VLBA band. 
We have also estimated individual errors cr^ for every source 
and frequency following two approache s: the one suggested by 
Fomalont ( 1999) and its modification bv lLee et all (^008). Then 
we compared the measured values of 6r with estimations of its 
error crgr for ev ery sub-band pai r. I f the o sr are estimated cor- 
rectly following lFomalon3 (1 1 9991) or lLee et al. (2008), we expect 
that for about 67% of th e sources os,- > Sr. We have found out 
that the lFomalontI (11999 ) method has provided meaningful error 
estimates for S, C, and X, while the L-bans errors were signif- 
icantly underestimated due to a very high signal-to- noise ratio 
of jet components in this band. At the same time, the lLee et"an 
('2008') method has significantly overestimated the errors in all 
the bands by 3 - 40 times. On the basis of these results, we made 
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Fig. 4. Distribution of the observed difference in distance be- 
tween the core and reference component measured at two sub- 
bands of X, C, S , and L bands. 



a decision to use a statistical method and estimate only a typical 
(Tr error for every band. An approach to adopting the same error 
for each target at a given frequency has an obvious weakness. 
However, we believe that this turns out to be the most reliable 
estimate of cr^ in view of the general problem of calculating the 
model-dependent error of a component position in VLBI. 

In the following, forX, C, and S bands we adopt a single typ- 
ical uncertainty of the distance between the core and reference 
component, which is the standard deviation of the distance dif- 
ference measurements made in the two sub-bands of each band, 
computed using all the sources in our experiment (Fig. |4). For 
the X band the adopted distance error is crx = 0.04 mas, for C 
band crc - 0.04 mas, and for S band os = 0.07 mas. The Ku 
band has not been divided into two independent sub-bands and 
for that band we adopt the same position uncertainty as for X 
band: ctk,, - (Tx - 0.04 mas. In the L band the mean difference 
between the distance measured at the lower sub-band (1 .4 GHz) 
and the higher sub-band (1.7 GHz) is -0.14 + 0.03 mas, which 
is significantly different from zero. We interpret this as a sign of 
the core shift effect, which should be more pronounced at lower 
frequencies, and use the uncertainty derived from the scatter of 
5 band measurements as an estimate of the L band error instead: 
(Tl-cts- 0.07 mas. We note that ctl may be underestimated. 

4. Discussion 

4. 1. Typically measured values of the shift 

In the 20 sources observed in our experiment we found the fol- 
lowing median values of the core shift with respect to its posi- 
tion at 15.4 GHz: 0.06 mas for 8.4 GHz, 0.08 mas for 8.1 GHz, 
0.22 mas for 5.0 GHz, 0.22 for 4.6 GHz, 0.76 mas for 2.4 GHz, 
0.72 mas for 2.3 GHz, 1.06 mas for 1.7 GHz, and 1.21 mas for 




Fig. 5. Distribution of the k parameter in the core position as 
a function of frequency fit: rdv) oc y^'^*. The mean value is 
k - 0.99 + 0.14, while the median is 0.82. One discrepant mea- 
surement on the histogram corresponds to the source 0952-)- 179 
where the value of k is not well constrained: k - 3.4 + 2.7. 



1.4 GHz measurements. See Table|4]for more details. It is diffi- 
cult to say how applicable these values are to the population of 
compact extragalactic radio sources as a whole. It is likely that 
the above values represent a higher extreme of the possible range 
of values, since the sources with a hi gh core shift were sel ected 
for our observations from the list of IKovalev et al.l (l2008l) . The 
median S/X band core shift value of 0.69 mas obtained for our 
sample in this study may be comp ared to the median value of 
0.44 mas from lKovalev et alj (|2008|) for a sample of 29 sources. 
These samples have 15 sources in common; see more detailed 
discussion in Section l43| 

O' Sullivan & Gabuzdaj (|200^ have measured the median 



core shift of 0.22 mas between 4.6 and 15.4 GHz using a small 
sample of three objects, none of which is present in our sample. 
This value is the same as the median 4.6/15.4 GHz shift in our 
sample (0.22 mas). A detailed investigation of the core shift ef- 
fect in a large complete sample of sources selected not on the 
basis of previous core shift measurements is needed to derive 
typical magnitude of the effect in a compact extragalactic radio 
source. 

Finally, we note that misalignment between the parsec-scale 
jet (resolved in our VLBI observations) and the subparsec scale 
jet (which we observe as the core) will reduce the frequency- 
dependent distance difference between the core and a reference 
feature in the parsec-scale jet. Therefore, the sources with the 
highest measured core shift may be those reflecting the true in- 
trinsic magnitude of the effect. 

4.2. Core position as a function of frequency 

According to iLobanovl d 19981) . the position of the core rdv) oc 
y-^/'' where the coefficient k - 1 if (i) the dominating absorp- 
tion mechanism is synchrotron self-absorption, the jet has 
a conical shape and ( Hi) is in equipartition. These assumptions 
must hold for the ultracompact jet region that we observe as the 
VLBI core and does not need to be correct all the way along the 
extended parsec-scale jet. 

To test the above assumptions, we fitted the observed dis- 
tances between the reference jet component and the core in each 
source with the function rdv) - a + /jv"'^* leaving the coeffi- 
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cients a, b, and k as free paramet ers. The Levenberg-Marquardt 
algorithm (e.g., Press et al. '2002*) was used to perform the non- 
linear least-squares fit. Figure |5] shows the distribution of the 
estimated k values. The mean value of the distribution is k - 
0.99 + 0.14. This suggests that, for a typical source in the sample 
in the frequency range of our observations, the above assump- 
tions about the structure of the inner jet (which we observe as the 
VLBI core) are generally correct. Specifically, we conclude that 
our observations are consistent with synchrotron self-absorption 
being the dominating opacity mechanism acting in parsec-scale 
core regions in the GHz frequency range. 

It should be noted that, since the core shift effect is most pro- 
nounced at lower frequencies, the best-fit values of k obtained 
heavily depend on the L-band measurements. At the same time, 
L-band data are also more likely to be affected by blending, 
an effect that could interfere with the core shift measurements 
dKov alev et alj |2008h . To test how sensitive our conclusions are 
to the L-band results, we have repeated the above fitting pro- 
cedure excluding all L-band data. The resulting values of k are 
on average 30 % lower than to those obtained with L-band data 
included. Naturally, a larger uncertainty is associated with the 
values of k obtained without L-band data. The mean values of 
the two distributions of best-fit k values are consistent at the 2cr 
level. 

The values of k presented in Fig. |5] are similar to 
those ob tained for 10 38 +528A 
1984 . iMarcaide et all (Il985h. 



fci_2 and ^2,3 is 



by Marcaide & Shapiro 



1998), 3C 309.1 by 'Ros & Lobano' 



3C345 
(l200l . 



3 (l2002h . 0850H-581 by 
iKovalev. Loba nov. & Pushkarev (200%J418-h546, 2007 -h777, 
BL Lac by [o ' Sullivan & Gabuzda (2009), and Mrk 501 by 
ICroke. O' Sullivan. & Gabu zda (2010). These values also agree 
with the results of I Yang et al. (2008), who estimates k from a 
statistical analysis of parsec-scale core sizes of ~ 3000 sources. 
The core size is expected to be directly rel ated to the observe d 
position of the core along the conical jet dUnwin et al.lll994b . 
The value of A: < 1 was obtained from analy sis of the core size i n 
3C 345 measured between 5 and 22 GHz bv lUnwin et"al](ll994l) . 
Values of k slightly more than unity were estimated for BL Lac 
by analyzing time lags between radio lightcurves obtained 
at diff'erent frequencies dBach et all [2006h . We also note that 
iLobanovl (Il998l) has found an indication that k changes with 
frequency up to v alues greater than 1 for the quasar 3C 309.1. 



by iLobanov 



iKadler etalj ("2004) found much higher values of k, indicating 
free-free absorption in the radio galaxy NGC 1052. The key 
difference between NGC 1052 and the sources observed in our 
sample is that this radio galaxy shows a two-sided parsec-scale 
jet. The inner part of the receding jet of NGC 1052 is lik ely to 
be obscured by a circumnuclear torus (IKadler et al .1120041) . 



4.3. Can the same function describe shifts at all frequencies 
between 1.4 and 15.4 GHz? 

To search for a possible change of k with frequency, following 
ILobanovl (1 1998h . we introduce a measure of core position offset 



Ar„ 



-\/k -\/k 
^1 - ^2 



between two frequencies (vi ,vi> v\ ), where Armas is the differ- 
ence in the apparent core positions at v\ and V2. Qrv can be used 
for assessing the possible variations of k. If two values Qrvi,2 
(measured between frequencies v\ and V2) and Q,v2,3 (between 
V2 and V3) are different, the relation between the corresponding 



Cl,2 



log 

log Q„ 



1,3 



Assuming ki^ — 1 for the lowest pair of frequencies (vi = 
1 .408 GHz and V2 = 1 .662 GHz for most cases), we derive the 
values of ^2,3, ^2,4, and so on. If the assumption about A: = 1 at 
L band is not correct, that would affect the resulting values of 
k\ however, the test will still be sensitive to relative changes in 
k between frequencies. The obtained values of k were consis- 
tent with A: = 1 at a 3cr level. Therefore we found no significant 
changes in k with frequency. 

The applicability of the « y^'/* relation to many ex- 
tragalactic sources has important consequences for the radio- 
optical reference frame alignment. In the particular case of 
k - 1, if the radio source is strongly core-dominated, its posi- 
tion measured using the group delay technique at any frequency 
(such as those currently used to d efine the the In t ernational 
Celestial Referen ce Frame ICRF , iMa et all Il998l iFev et al.1 
,2004. and ICRF2, iFev et alJlIoToh corresponds t o the jet base , 
not to the radio core position at a given frequency (iPorcas l2009h . 
This is due to the additional time delay introduced by the core 
shift across the observing band (Porcas 2009). The optical emis- 
sion of an AGN jet is expected to originate much closer to the 
jet base then the radio core position. Therefore, if the structure 
of the radio source is dominated by the core, there is no need 
to introduce core shift corrections to the radio source positions 
determined from group delays before matching them with op- 
tical positions. This is not true, however, for radio so urce posi- 
tions determined from phase referencing ( Porcas! E00 9). If a ra- 
dio source has distinct jet components comparable in brightness 
to the core (such as many sources considered here; Fig. |T}, the 
group delay would not be a good estimator of the jet base posi- 
tion. The source structure needs to be taken into account prior to 
comparing the radio and optical source positions in such cases. 

4.4. Fitting results assuming opacity due to synchrotron 
self-absorption 

As discussed above, our observations are consistent with the 
assumption of synchrotron self-absorption dominated opacity 
across the frequency range of this VLBI experiment. We re- 
peated the fit of the function rdv) - a + Zjv^'^* to the observa- 
tional data, fixing the value ofk- 1, which is our best estimate. 
The results are presented in Table [Sjand Fig. |3] The values of a 
and b are given in TablejS] The projected core distance in parsecs 
is 



Da , 

rdv) [pc] = a [mas] + b [mas • GHz] • 



lad 



1 -^z 



(GHz) ] 



where a [mas] and b [mas ■ GHz] are the coefficients obtained 
from observations. Da is the angular size distance to the source 
in parsecs, A^^d - 206264800 is the number of milliarcseconds 
in one radian, z the source redshift, and Vem (GHz) the emission 
frequency in the source frame. 

We have investigated deviations of the measured distances 
from the best-fit curve presented in Fig. |3]for each frequency. 
In all cases the mean difference between the measurement and 
the model curve is consistent with zero, and its standard devia- 
tion is consistent within a factor of two with the adopted typical 
measurement uncertainties. 

It can be seen from Fig.|3]that, while the rdv) cc y"' law gen- 
erally provides a good fit to the data over the whole observed fre- 
quency range with no systematic deviations, there are indications 
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Table 5. Distance of the reference jet component from the core as a function of frequency: fit results (r - a + blv) 



Source 


a (mas) 


b (mas ■ GHz) 


a (pc) 


b (pc ■ 


GHz) 


P.A.° 


0148+274 


10.32 ± 0.05 


-2.48 ±0.19 


86.75 


±0.42 


-47.11 


±3.61 


-40.9 


0342+147 


7.15+0.03 


-1.08 ±0.10 


61.14 


±0.26 


-23.61 


±2.19 


-88.5 


0425+048 


18.58 + 0.06 


-2.20 ±0.19 


115.05 


±0.37 


-20.66 


± 1.78 


-100.7 


0507+179 


4.52 ± 0.05 


-1.73 ±0.19 


24.74 


±0.27 


-13.41 


± 1.47 


-100.3 


0610+260 


6.34 ± 0.06 


-3.69 ±0.20 


41.62 


±0.39 


-38.28 


±2.07 


-94.5 


0839+187 


12.21 +0.09 


-3.67 ±0.31 


102.76 


±0.76 


-70.18 


±5.93 


15.2 


0952+179 


14.46 + 0.06 


-1.60 ±0.19 


123.39 


±0.51 


-33.83 


±4.02 


-3.5 


1004+141 


11.25 ±0.04 


-2.53 ±0.13 


90.53 


±0.32 


-75.47 


±3.88 


132.2 


1011+250 


6.89 ± 0.05 


-2.10 ±0.15 


58.96 


±0.43 


-47.37 


±3.38 


-107.3 


1049+215 


8.33+0.02 


-1.82 ±0.09 


70.30 


±0.17 


-35.33 


± 1.75 


108.0 


1219+285 


6.95 + 0.04 


-2.38 ±0.16 


19.06 


±0.11 


-7.58 


±0.51 


109.8 


1406-076 


7.25 ±0.01 


-1.25 ±0.05 


61.90 


±0.09 


-26.61 


± 1.06 


-103.5 


1458+718 


24.57 ± 0.05 


-2.38 ±0.18 


192.21 


±0.39 


-35.45 


±2.68 


163.7 


1642+690 


10.62 ±0.05 


-1.68 ±0.19 


78.01 


±0.37 


-21.61 


±2.44 


-167.0 


1655+077 


8.50 ±0.02 


-1.37 ±0.08 


57.64 


±0.14 


-15.06 


±0.88 


-42.9 


1803+784 


6.83 ± 0.04 


-1.04 ±0.14 


48.21 


±0.28 


-12.33 


± 1.66 


-96.0 


1830+285 


6.40 ± 0.06 


-3.15 ±0.24 


42.50 


±0.40 


-33.34 


±2.54 


-38.5 


1845+797 


7.23 ± 0.09 


-2.29 ±0.35 


7.76 


±0.10 


-2.59 


±0.40 


-38.0 


2201+315 


6.78 ± 0.03 


-3.22 ±0.15 


29.83 


±0.13 


-18.39 


±0.86 


-135.3 


2320+506 


8.52 ±0.03 


-1.15 ±0.09 


71.76 


±0.25 


-22.07 


± 1.73 


-135.8 



Column designation: Col. 1 - lAU source name (B1950), Col. 2 and 3 - coefficients of the best-fit curve r - a + b/v with their uncertainties. 
Col. 3 and 4 - same coefficients converted to the projected linear scale. Col. 5 ~ median (across all frequencies) position angle of the reference jet 
component with respect to the core; it marks the direction along which the core shift was measured. 





2 


4 - 




2 


2 - 


o 


2 


- 


CO 






shift 


1 


8 - 


Core 


1 


6 - 




1 


4 - 




1 


2 - 




1 


- 




0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 
Core flux density ratio 

Fig. 6. Ratio of the X-S ban d core position shift measured at 
2002 bv ,Kovalev et alJ (12008 ) to the one measured in 2007 (this 
paper) as a function of the X band core flux density ratio at these 
two epochs. The core shift ratio is defined to be greater than 
unity, i.e., for each pair of measurements, the higher value of the 
shift is in the numerator of the ratio. The one source (not shown 
in this plot), which exhibited a large core shift change without 
a major change in flux density is W Com (1219+285); it has a 
core shift ratio of 3.67 + 0.36 and core flux ratio of 1.19 + 0.08. 



of possible local deviations. This may be due to local deviations 
from equipartition or conical shape of the jet. Another possibil- 
ity is that the scatter in distance measurements at two sub-bands 
used for error estimation does not reflect the total measurement 



uncertainty. This may happen if there is a systematic factor that 
influences distance measurements at both sub-bands in a similar 
way. However, after a careful review of our analysis procedures 
we could not identify any such factor. 

4.5. Search for core shift variability 

The new measurements of core shifts between 2.3 and 8.4 GHz 
may be compared to the previous measurements at similar fre- 
quencies obtained using the RDV data bv lKovalevet alJ ( l2008h . 
As expected for a big enough sample of variable sources, the 
mean difference between the measurements obtained at two 
epochs (2002 for RDV global VLB! data and 2007 for this 
VLBA dataset) is consistent with zero: 0.01 + 0.08 mas. The 
mean absolute value of this difference is 0.22 + 0.05 mas, which 
is a measure of a typical long-term core-shift variability in the 
studied sample within the accuracy of our measurements. 

According to formula (1) in Kova lev et al.l ([2008), if jet ori- 
entation and velocity remain unchanged, an increasing value of 
core shift (due to an increase in particle density and/or jet mag- 
netic field strength) should be associated with an increase in the 
flux density of the core. Figure|6]shows the core shift ratio versus 
X band core flux density ratio at two epochs (2002 and 2007). No 
significant correlation is found between the core shift and flux 
ratios at different epochs. This may reflect that there is no sin- 
gle dominant physical reason for core flares and subsequent core 
shift variations, despite the suggestion by Kovalev e t al. (200^. 
A higher accuracy of core-shift ratio measurements is required 
to probe this in detail. 

5. Summary 

Dedicated VLBA observations of twenty AGN jets showing 
large frequency-dependent core shifts were conducted at nine 
frequencies in the 1.4-15.4 GHz range. Significant core shifts 
have been detected and confirmed in all observed sources. The 
effect is more pronounced at lower frequencies. The median 
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value of the core shift in the observed sources is 1.21 mas if 
measured between 1.4 and 15.4 GHz and 0.24 mas between 5.0 
and 15.4 GHz. The core position, r^, shift as a function of fre- 
quency, V, is consistent with the rdv) oc y^' pattern expected 
for a purely synch rotron self-absorbed conical jet in equip arti- 
tion dBlandford & Konigl 1979, Lobanov 1998). The mean value 
of the coefficient k in the rdv) oc y^'^* relation was found to 
he k = 0.99 + 0.14, while the median is 0.82. These results 
support the interpretation of the parsec-scale core as a contin- 
uous Blandford-Konigl type jet with smooth gradients of physi- 
cal properties (including opacity) along it. No systematic change 
with frequency of the power law index in the rdv) relation has 
been convincingly detected. However, some local changes might 
be present in a few sources, especially at higher frequencies. 

The general applicability of the rdv) oc y^' relation to all 
the observed sources is a promising indication that, if a radio 
source is strongly core-dominated, its positions obtained using 
group delays (such as those currently used to define the 1CRF2) 
may be compared to optical counterpart positions direct ly, with 
no ne ed to apply a correction for the source core shift (iPorcasI 
I2OO9I) . This may be important in the era of the future space- 
based optical astrometry missions aiming for //as-level accuracy. 
A correction for the core shift in the reference source(s) is still 
required if a high-precision absolute position of a celestial ra- 
dio source (or tracked spacecraft) is to be determined through a 
phase-referencing VLBl experiment. 

While in this paper we have tried to concentrate on the ob- 
servational results, independent of specific assumptions about 
jet geometry and Lorentz factor, in future work based on the 
dataset presented here we plan to estimate geometry, magnetic 
field strength, and total (kinetic plus magnetic field) power of 
the flows and to relate the observed shifts to the properties 
of the central black hole and broad-line region in its vicinity. 
Polarization information will be used to constrain the physical 
interpretation of core shift measurements and to investigate the 
effect of the core shift on Faraday rotation measurements. 

It would be important to confirm and investigate possible 
changes in the coefficient k in the rdv) 0^ y^'/* relation with fre- 
quency, which may be hinted at in our observations. To achieve 
greater accuracy in measuring k one would need to obtain more 
independent core position measurements across the whole fre- 
quency range accessible to VLBl. The use of phase-referenced 
observations may also be helpful for minimizing ambiguities in- 
evitably associated with the model fitting procedure and quanti- 
fying position-measurement errors more accurately. 

Calibrated Mv-data and models described in this paper are 
available from the first author (KVS) by request. 
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Fig. 1. Naturally weighted CLEAN images of the observed sources between 1.4 and 15 GHz. The lowest contour value 'clev' is 
chosen at four times the rms noise, the peak brightness is given by 'max' . The contour levels increase by factors of two. The dashed 
contours indicate negative brightness. The beam's full width at half maximum (FWHM) is shown in the bottom left comer of the 
images in gray. An epoch of observation is shown in the bottom right corner. Red and blue spots indicate the positions and sizes 
(FWHM) of Gaussian model components for the core and the jet features, respectively. 
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